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We evaluate f n /m p , fx/i^p, and vn^j '(f$m p ), extrapolated to physical 

quark mass, zero lattice spacing and infinite volume, for lattice QCD with Wilson 
quarks in the valence (quenched) approximation. The predicted m^j '(f<f,m p ) differs 
from experiment by less than its statistical uncertainty of approximately 15%. The 
other three constants are 10% to 20% below experiment, equivalent to between one 
and two times the corresponding statistical uncertainties. 
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In a recent paper |IJ we presented lattice QCD predictions for the masses of eight low- 
lying hardrons, extrapolated to physical quark mass, zero lattice spacing, and infinite vol- 
ume, using Wilson quarks in the valence (quenched) approximation. The masses we found 
were within 6% of experiment, and all differences between prediction and experiment were 
consistent with the calculation's statistical uncertainty. We argued that this result could be 
interpreted as quantitative confirmation of the low-lying mass predictions both of QCD and 
of the valence approximation. It appeared to us unlikely that eight different valence approx- 
imation masses would agree with experiment yet differ significantly from QCD's predictions 
including the full effect of quark- ant iquark vacuum polarization. 

We have now evaluated the infinite volume, zero lattice spacing, physical quark mass limit 
of fn/m p , fx/rrip, l/f p , and ra^j '(f^m p ). To our knowledge there have been no previous 
systematic calculations of this physical limit of lattice meson decay constants. A review 
of earlier work is given in Ref. 0. The predicted m^/^fymp) differs from experiment by 
less than its statistical uncertainty of approximately 15%. The other three constants are 
10% to 20% below experiment, equivalent to between one and two times the corresponding 
statistical uncertainties. 

The decay constant predictions we obtain may appear to be marginally in conflict with 
our earlier mass results. A simple argument suggests, however, that valence approximation 
decay constants actually should lie at least some amount below the predictions of the full 
theory and therefore, presumably, below experiment 0. The valence approximation may be 
viewed as replacing the momentum and frequency dependent color dielectric constant arising 
from quark-antiquark vacuum polarization with it zero-momentum limit |4j]. The valence 
approximation might thus be expected to be fairly reliable for low-lying baryon and meson 
masses, which are determined largely by the low momentum behavior of the chromoelectric 
field. The effective quark charge used in the valence approximation may be thought of as 
the product of the zero-momentum dielectric constant and the quark's true charge. The 
valence approximation's effective quark charge at higher momenta can then be obtained 
from the low momentum charge by the Callan-Symanzik equation. The leading coefficient 
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(3q entering the Callan-Symanzik equation in the valence approximation is larger than (5q of 
the full theory, which is reduced by dynamical quark vacuum polarization. It follows that 
the quark charge in the valence approximation will fall faster with momentum than it does 
in the full theory At sufficiently short distance the attractive quark- ant iquark potential 
in the valence approximation will therefore be weaker than in the full theory. Meson wave 
functions in the valence approximation will then be pulled into the origin less than in the full 
theory and will be smaller at the origin. Since meson decay constants are proportional to 
quark- ant iquark wave functions at the origin, decay constants in the valence approximation 
will be smaller than in the full theory. 

Thus while Ref. [IJ confirms the valence approximation for the low momentum behavior 
of the chromoelectic field governing the binding energy of low-lying hadrons, the present cal- 
culation shows that the valence approximation appears to become less reliable, as expected, 
for the higher momenta contributing to meson wave functions at the origin. Our result 
suggests that valence approximation predictions for the hadronic matrix elements entering 
weak decays will also tend to be somewhat smaller than their values in the full theory. 

The calculations described here were done on the GF11 parallel computer at IBM Re- 
search |J and use the same collection of gauge configurations and quark propagators gener- 
ated for the mass calculations of Ref. Q . The full set of mass and decay constant calculations 
took approximately one year to complete. GF11 was used in configurations ranging from 
384 to 480 processors, with sustained speeds ranging from 5 Gflops to 7 Gflops. With the 
present set of improved algorithms and 480 processors, these calculations could be repeated 
in less than four months. 

The normalization we adopt for pseudoscalar and vector decay constants in continuum 
QCD is < 0| X- M |P(p, j) > = ip^f) and < 0| Jj\V(p, e, j) > = e^rrijFj, for states normalized 
by <p\q> = (2tt) 3 p°5(p — q), where Jj 1 and Jj^ are vector and axial vector flavor-SU(3) 
currents and j is an adjoint-representation flavor index running from 1 to 8. Assuming exact 
isospin symmetry, we have fi = f^, for i = 1, ...3, fi = fx, for i = 4, ...7, and Fi = m p / f p , 
for i = 1,...3. In the valence approximation we have, in addition, F 8 = 3m^/(\/2/^)- 
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For simplicity, we will also use the names F p and for F\ and F$, respectively. Our 
normalization gives f n the value (93.15 ± 0.11)MeV. 

The hadron mass calculation of Ref. was done using gaussian smeared quark and 
antiquark fields defined in Coulomb gauge. Smeared fields have, therefore, also been adopted 
for the present calculation. The smeared field 4> r (x, t) is J2yG r (x — yjtpiy, t) for a smearing 
function G r (z) of (y / 7rr)~ 3 exp( — \z\ 2 /r 2 ), and (p r (x, t) is defined correspondingly. We take the 
smeared fields (f>o(x) and o (x) to be ip{x) and ip(x), respectively. From these fields, define 
the composite field Jj r and smeared lattice currents Jj r and J_y to be r 7 5 Aj0 r , r 7 At Aj^ r 
and <j) r , y 5 'y fl \j<j) r , respectively, where the Xj are an orthonormal basis for the flavor-SU(3) 
Lie algebra with normalization tr(XjXj) = 1/2. 

Define the constants ZU* r , Z^F r , and Zj r Y r by the requirement that for a lattice with 
time direction period L and a time t, L ^> t ^> 1, the correlation function J2s < 
[Jj r ,(x, t)}^ Jj r (0, 0) > approaches Z^ r ex"p[—m^t)^ J2x < [Jjr'(%i t)]*Jjr(®i 0) > approaches 
Z£F r exp(— mft), and J2s < [^j r '(^^)] t ^jr(0; 0) > approaches ZjP£e:rp(— mjt). Here m^ 5 and 
are pseudoscalar and vector masses, respectively. 

Measured in units of the lattice spacing a with any choice of the smearing radius 
r, the decay constant fja then is given by QhimjaZj^, )~^zfZ^£ , and Fja is given by 
2a (rrijaZy^)~2 zj Zj^ . The coefficients zf and zj 7 are finite renormalizations chosen so 
that the lattice currents z^J^ and zj J^ r approach the continuum currents j| M and Jj 4 , 
respectively, as the lattice spacing approaches zero. A mean-field theory improved pertur- 
bation expansion || gives zf the value (1 — 3kj/4k c )[l — 0. 310^(1/(2)] and zj the value 
(1 — 3kj/4k c )[l — 0.82ams(l/a)], where kj is the hopping constant corresponding to the mass 
of the quark and antiquark for a meson with flavor j, assumed here to have m q = m-g, and k c 
is the critical hopping constant at which the pion's mass becomes zero. The decay constant 
for a meson with m q ^ m-g will be discussed below. The renormalization constants zf and 
zj are often given the "naive" value 2k j. Even to leading order in perturbation theory, 
however, a better choice turns out to be (1 — 3kj/4k c ) 0. 

Table | lists the lattice sizes, parameter values, sweeps skipped between gauge config- 
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urations, and number of configurations used in the ensembles from which decay constants 
were calculated. Gauge configurations were generated with the Cabbibo-Marinari-Okawa 
algorithm. A discussion of the algorithms by which quark propagators were found is given 
in Ref . |I| . Hadron masses were determined by fits to hadron propagators constructed from 
the quark propagators. The coefficients > ^22; ^j02^ anc ^ were then extracted 
from fits to current-current and pseudoscalar-current correlation functions, and combined 
with both perturbative and naive renormalization to find two different sets of values for 
f n and F p . Statistical errors on these quantities and all other parameters which we have 
determined were found by the bootstrap method. A more detailed discussion of our fits and 
error analysis will be given elsewhere [0]. 

Percentage changes in decay constants going from 16 3 x 32 to 24 3 x 32, at (5 of 5.7, are given 
in Table |I[ These changes are the same for both perturbative and naive renormalization. 
All of the differences appear to be of marginal statistical significance and may therefore 
best be viewed as upper bounds on the volume dependence of our results. It appears quite 
likely that for the range of k, (3, and lattice volume we have examined the errors in valence 
approximation decay constants due to calculation in a finite volume L 3 are bounded by an 
expression of the form Ce~ L ^ R , with a coefficient R of the order of the radius of a hadron's 
wave function. At f3 of 5.7 for the k we considered, R is thus typically 3 lattice units or less. 
We therefore expect that the differences between decay constants on a 16 3 volume and those 
on a 24 3 volume are nearly equal to the differences between 16 3 and true infinite volume 
limiting values. 

At the largest k on each lattice, the ratio m^/nip is significantly larger than its experi- 
mentally observed value of 0.179. Thus to produce decay constants for hadrons containing 
only light quarks, our data has to be extrapolated to larger k or, equivalently, to smaller 
quark mass. We did not calculate directly at larger k both because the algorithms we used 
to find quark propagators became too slow and because the statistical errors we found in 
trial calculations became too large. 

To extrapolate F p and f n to small quark mass, we first determined the critical hopping 
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constant k c , for each lattice and /3, at which m n becomes zero |T|. Defining the quark mass 
in lattice units m q a to be l/(2k) — l/(2k c ), we found f n a and F p a, both for perturbative 
renormalization and for naive renormalization, to be nearly linear functions of m q a over the 
entire range of k considered on each lattice. Figure [I] shows f n and F p , given by perturbative 
renormalization, as functions of m q . Data is shown from all lattices of Table [I] except 
24 3 x 32 at (3 = 5.7. For convenience, we show all hadron masses in units of the physical rho 
mass, m p (m n ), given by m p evaluated at the "normal" quark mass m n which produces the 
physical value of m n /m p . The quark mass m q in Figure [l] is shown in units of the strange 
quark mass m s . The value of m s for each lattice and (3 is found in Ref. Jl| by requiring 
"v[( m n + m s)/2]/m p (m n ) to be equal to the physical value of m K /m p . The straight lines 
are fits to the m q dependence of f n and F p at the three smallest quark masses in each 
data set at fixed (3. The fits in Figure [I] appear to be quite good and provide, we believe, 
a reliable method for extrapolating decay constants down to light quark masses. With 
naive renormalization, /„. and F p fit straight lines in m q about as well as the perturbatively 
renormalized data of Figure [I]. 

The linear fits of Figure [TJ permit the determination of fx and in addition to f n and 
F p . For a pion composed of a quark and antiquark with mass m q ^ m q , Figure [I] suggests 
fn = (y q m q + a q m q + (3. Charge conjugation invariance then gives a q = a q . It follows that the 
kaon, which is a pion with, say, m q = m s and m q = m n , will have the same decay constant 
as a pion composed of a single type of quark and antiquark with m q = m q = (m s + m n )/2. 
On the other hand, the linear fits of Figure [l] permit F p to be extrapolated to the point 
m q = m q = m s which, in the valence approximation, gives F^. 

The ratios f n /m p , /^/m p , F p /m p and F^/rrip for physical quark masses we then extrapo- 
lated to zero lattice spacing. The physical lattice volume was held nearly fixed as the lattice 
spacing was taken to zero. For Wilson fermions the leading lattice spacing dependence in 
mass ratios is expected to be linear in a. Figure ^| shows decays constants with perturbative 
renormalization along with linear fits to m p a. The quantity m p a may be viewed as the lattice 
spacing a measured in units of the physical rho Compton wavelength, l/m p . The vertical 
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bars at m p a of 0, offset slightly for visibility, are the extrapolated predictions' uncertainties 
determined by the bootstrap method, and the horizontal bars are corresponding experimen- 
tally observed values. The calculated points in Figure ^ are for the lattices 16 3 x 32, 24 3 x 36 
and 30 x 32 2 x 40. The values of (3 for these lattices were chosen so that the physical volume 
in each case is nearly the same. For lattice period L, the quantity m p L is respectively, 9.08 ± 
0.13, 9.24 ± 0.19 and, averaged over three directions, 8.67 ± 0.12 (TJ. Extrapolations to zero 
lattice spacing similar to those shown in Figure § were also done for naive renormalization. 

The continuum ratios we found in finite volume were then corrected to infinite volume 
by an adaptation of the method used in Ref. [0] to correct finite volume continuum mass 
ratios to infinite volume. The difference between a unitless decay ratio found on the lattice 
16 3 x 32, at (3 of 5.7, and the same ratio found on the lattice 24 3 x 32, at f3 of 5.7, we took as 
a finite lattice spacing approximation to the difference between the continuum decay ratio 
in a box with period having m p L of 9 and the continuum decay ratio in infinite volume. The 
error in this procedure, for perturbatively renormalized decay constants, can be estimated 
to be 4 % or less as follows. For perturbative renormalization, the value of each decay 
ratio on the lattice 16 3 x 32 at (3 of 5.7 differs from its continuum limit by at most 75 % 
of the continuum limit. Moreover, as we argued earlier, the change in each decay ratio, at 
(3 of 5.7, between 16 3 x 32 and 24 3 x 32 should be nearly the same as the corresponding 
change between 16 3 x 32 and infinite volume. Combining these two pieces of information, 
we expect that with a relative error of 75 % or less, the change in any decay ratio between 
16 3 x 32 and 24 3 x 32 at (3 of 5.7 should be the same as the change between the continuum 
decay ratio in a box with period having m p L of 9 and the corresponding continuum ratio 
in infinite volume. Since the changes we found in decay ratios, for physical quark mass, 
between 16 3 x 32 and 24 3 x 32 are all less than 5 %, the error in using these differences 
as estimates of corresponding continuum differences between m p L of 9 and infinite volume 
should of the order of 75 % of 5%, which is 3.75 %. The corresponding uncertainty for naive 
renormalization is 10.0 % due to a larger change in these values between 16 3 x 32 at (3 of 
5.7 and the continuum limit. 
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The ratios f n /m p , fx/rrip, F p /m p and F^/rrip, for both naive and perturbative renor- 
malization, extrapolated to zero lattice spacing with m p L fixed at 9, and then corrected to 
infinite volume are shown in Table |T1|. The errors shown for infinite volume ratios are statis- 
tical only and do not include the estimates we have just given for the systematic error in our 
procedure for making infinite volume corrections. Most of the perturbatively renormalized 
results lie below the corresponding ratios with naive renormalization. For all ratios except 
F<p/m p in finite volume, however, the two different renormalizations give results which are 
statistically consistant. Where there is a disagreement, the perturbatively renormalized ra- 
tios are more reliable estimates of the true infinite volume continuum limit of the valence 
approximation. For extrapolations done from lattice spacings which are sufficiently small, 
the two different renormalization procedures should give the same limiting results. The 
agreement between the two different methods for three of four parameters tends to support 
the reliability of our extrapolations to the continuum limit. 

The predicted infinite volume ratios in Table [TT1] are statistically consistent with the finite 
volume ratios. The main consequence of the correction to infinite volume is an increase in 
the size of the statistical uncertainty in each prediction. 

Three of the four finite volume volume continuum ratios obtained by perturbative renor- 
malization lie below experiment by between one and two standard deviations. The larger 
physical volume over which infinite volume meson wave functions can spread should give 
infinite volume decay constants some amount smaller than those calculated in finite vol- 
ume. Thus the true disagreement between infinite volume f n /m p , f K /m p , and F p /m p and 
experiment should be at least as great than the disagreement we find for the finite volume 
predictions. The statistical significance of the disagreement between our inifinite volume 
results themselves and experiment, however, is weaker than for the finite volume predictions 
as a result of the large size of the error bars on the infinite volume predictions. The larger 
error bars of the naive renormalized decay constants also weakens the significance of their 
disagreement with experiment. 
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FIGURES 

FIG. 1. The decay constants F p and /„-, with perturbative renormalization, in units of the 
physical rho mass m p (m n ), as functions of the quark mass m q , in units of the strange quark mass 
m s . 

FIG. 2. Perturbatively renormalized decay constants as functions of the lattice spacing a, in 
units of l/m p . The error bars near zero lattice spacing are uncertainties in the extrapolated ratios, 
and the horizontal lines represent experimentally observed values. 
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TABLES 



lattice 




k 


skip 


count 


16 3 x 32 
24 3 x 32 
24 3 x 36 
30 x 32 2 x 40 


5.7 
5.7 
5.93 
6.17 


0.1600 - 0.1675 
0.1600 - 0.1675 
0.1543 - 0.1581 
0.1500 - 0.1532 


2000 
4000 
4000 
6000 


219 
58 
217 
219 


TABLE I. Configurations analyzed. 


decay 




k 




change 


u 




0.1600 
0.1650 
0.1663 
0.1675 




3.510% 
-l.lli;i% 

-3.i±i:J% 

-6.2l^% 


fp 




0.1600 
0.1650 
0.1663 
0.1675 




7-lt^% 

4.9l?;°% 

4.010% 
-2.3l^|% 



TABLE II. Changes in perturbatively renormalized decay constants from 16 3 x 32 to 24 3 x 32 
at = 5.7 
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decay 


renorm. 


finite volume 


infinite volume 


obs. 


f / Tn 
JTT 1 "<y> 


UCl U LI J. U. 

naive 


U.lUZ_g Q 2 g 

n nn 7 +0.030 
u - uy ' -0.032 


091+^ 
u.u»i_q q 36 

o-o82±8:gg 


0.121 


f js hrn 
J K / "<-/> 


YJ \Zi Ij LLJ. \ J . 

naive 


119+^1 

U.11»_Q Q 15 

U.lOU_ 022 


u.iiu_g Q24 

n i 97 +0.030 
u - iz ' -0.032 


0.148 


F p /m p 


perturb, 
naive 


77+0 .017 
U- 1 ' ' -0.020 

n -, ql +0.033 

u.iyi_ 0031 


n i 70+O.O39 
U.l(O_ 0031 

n 9fl1 +0.035 
U./Ul_o 059 


0.199 


F^/nrLp 


perturb, 
naive 


u - zi ' -0.022 

2 70+ - 021 

u.z/u_ 029 


u.zoo_ 051 


0.219 



TABLE III. Calculated values of meson decay constants extrapolated to zero lattice spacing in 
finite volume, then corrected to infinite volume, compared with observed values. 
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